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The hyperfine structure of the interstitial muonium (Mu) in rutile (TiC> 2 , weakly n-type) has been 
identified by means of a muon spin rotation technique. The angle-resolved hyperfine parameters 
exhibit a tetragonal anisotropy within the ab plane and axial anisotropy with respect to the (001) 

(c) axis. This strongly suggests that the Mu is bound to O (forming an OH bond) at an off- 
center site within a channel along the c axis, while the unpaired Mu electron is localized around the 
neighboring Ti site. The hyperfine parameters are quantitatively explained by a model that considers 
spin polarization of the unpaired electron at both the Ti and O sites, providing evidence for the 
formation of Mu as a Ti-O-Mu complex paramagnetic state. The disappearance of the Mu signal 
above ~10 K suggests that the energy necessary for the promotion of the unpaired electron to the 
conduction band by thermal activation is of the order of 10 1 meV. These observations suggest that, 
while the electronic structure of Mu (and hence H) differs considerably from that of the conventional 
shallow level donor described by the effective mass model, Mu supplies a loosely bound electron, 
and thus, serves as a donor in rutile. 

PACS numbers: 61.72.Hh, 71.55.Ht, 76.75.+i 


I. INTRODUCTION 

Hydrogen (H) is a ubiquitous impurity in a wide vari¬ 
ety of semiconductors, including elemental (e.g., Si) and 
binary compound (e.g., GaAs) materials, which comprise 
the primary basis for current industrial applications. In 
these covalent systems, H is known to be amphoteric, 
forming an acceptor or donor level in n- or p-type ma¬ 
terials, respectively^ - — . Meanwhile, a recent theoretical 
prediction that H could behave as an independent elec¬ 
tron donor to induce n-type conductivity in ZnO^, along 
with subsequent support from various experiments in¬ 
volving muon spin rotation (/iSR)&£, electron param¬ 
agnetic resonance (EPR), electron-nuclear double reso¬ 
nance (ENDOR) measurements^, and additional conven¬ 
tional techniques^ - —, have prompted the development of 
a generalized hypothesis of H as a shallow donor dopant 
in wide gap semiconductors^ - — . 

As in the case of ZnO, /iSR studies have made a sig¬ 
nificant contribution to this subject by providing spec¬ 
troscopic information on muonium (Mu; an analog of a 
neutral H atom where the proton is replaced by a pos¬ 
itive muon) observed in a wide variety of semiconduc¬ 
tors. While the dynamical aspects (e.g., diffusion) of 
Mu and H may differ considerably, because of the rel¬ 
atively low mass of Mu (m p ~ m p / 9, with m p being 
the proton mass), the local electronic structure of Mu is 
virtually equivalent to that of H if a small correction to 
account for the difference in the reduced mass (^0.5%) 
is made. The recent discovery of novel Mu states hav¬ 
ing extremely small hyperfine parameters (^10 -4 times 
smaller than that of Mu in vacuum) and low ionization 
energy (^lO 1 meV) in several compound semiconductors, 
including CdS^, ZnSeii, InNi^, and Gahti^, also supports 


the hypothesis that Mu (and hence H) can act as a donor 
in these compounds. 

The generalized hypothesis also predicts shallow donor 
H in various metal oxides, most notably in rutile Ti02, 
where the origin of the unintentional n-type conductivity 
exhibited by the as-grown crystal is the focus of intense 
research. While the electrical and optical properties of 
Ti 02 are subject to both intrinsic and extrinsic defects, 
this material has a rich variety of potential applications in 
both electronic and opto-electronic devises, and detailed 
understanding of the H behavior is crucially important 
to the material functionality in this context. 

As shown in Fig. QJa), rutile has a tetragonal struc¬ 
ture with Ti 4+ ions occupying the body-center position, 
where each Ti ion is coordinated by six ligand oxygen 
atoms (O 2- ) comprising a slightly distorted octahedron. 
The TiOe octahedra are linked by corner-sharing to form 
O chains parallel to the c axis. Regions of low electron 
density exist between the chains, forming the so-called 
“c channel”. Earlier studies have shown that interstitial 
H resides in the c channel, leading to the formation of 
an OH bond (with a length of 0.109 nm) perpendicu¬ 
lar to the c axisii. The open channels cause anisotropy 
in the diffusion process, and may allow fast diffusion of 
smaller ions parallel to the tetragonal axis. It has re¬ 
cently been reported that the H diffusion can be strongly 
enhanced via irradiation of the crystal by infrared light 
that matches the OH bond stretch modei^. 

Here, we report on the electronic structure of Mu in 
single-crystalline rutile. Although the existence of a Mu 
state with an extremely small hyperfine (HF) parameter 
has been reported elsewher o 19 i 2Q , little is known about 
the local electronic structure. In this study, we show 
that the HF parameters exhibit strong anisotropy within 
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FIG. 1: (Color online) (a) Schematic illustration of rutile 
Ti 02 crystal structure, showing a unit cell with a — b — 
0.4593 nm and c = 0.2959 nm. (b) Atomic position of in¬ 
terstitial hydrogen (H), where H forms OH bonds (0.109 nm) 
at four inequivalent H sites, (c)-(e) Signs of hyperhne pa¬ 
rameters ( Ai ) expected for the magnetic dipolar field (Rdip) 
generated by an electron spin localized on Ti atom under an 
external field ( B ). (see text). 


the ab plane as well as along the c axis. The angular de¬ 
pendence of the HF parameters is perfectly in line with 
that expected for Mu occupying the sites that are com¬ 
mon to the interstitial H and forming OMu bonds. More 
interestingly, the relative signs of these parameters indi¬ 
cate that the HF interaction is predominantly determined 
by the magnetic dipolar fields generated by the unpaired 
Mu electron located away from the muon site. This is 
qualitatively consistent with a recent report on the H- 
related paramagnetic center, in which the HF interaction 
is attributed to the Ti 3+ [OH] “ complex state with the d 
electron being located on the Ti ion nearest to the OH 
basest Meanwhile, a detailed comparison of the Mu and 
H HF structures reveals a distinct difference between the 
complex centers, where a contribution from the residual 
spin polarization at the nearest neighboring O site (an¬ 
tiparallel to that of the Ti ion) is suggested in the Mu 
case. Through consideration of the temperature (T) de¬ 
pendence of the Mu yield, we also show that the unpaired 
electron bound to the Mu complex state requires small 
energy for promotion to the conduction band. Based on 
these observations, we argue that interstitial Mu/H is one 
of the primary origins of unintentional n-type conductiv¬ 
ity in rutile TiCU- This is despite its complex electronic 
structure, which differs significantly from that for the 
conventional effective-mass-like donor. 


II. EXPERIMENT 

A conventional fi SR experiment was conducted at the 
former Booster Meson Facility in Tsukuba (at KEK) and 


at J-PARC MUSE. A 4-MeV muon beam that was almost 
100 % spin-polarized parallel to the beam direction was 
implanted into single-crystalline rutile wafers [30 mm x 
30 mm, 1 . 0 -mm thick, grown using the Verneuil process], 
which were “as received” from a local vendor (Furuuchi 
Chemical Co.). The definitions of the azimuthal and po¬ 
lar angles (</>, 0 , respectively) for the external magnetic 
field ( B ) are shown in Fig. Oja). For the ^-dependence 
measurements, a wafer having a normal axis orientation 
of ( 001 ) (|| c) was used, where the initial muon spin po¬ 
larization (P M ) was parallel to c, and the a axis was ro¬ 
tated contrary to the B direction within the ab plane (i.e., 
0 = 90°). Meanwhile, the 0-dependence was measured 
by tilting the B direction towards the c axis (except for 
0 = 0 , for which another ( 100 ) crystal was used), where 
the angle was tuned by an additional field parallel to the 
c axis. In both cases, the P-parallel component of the HF 
parameter was observed in terms of the frequency shift 
of the satellite signals. 


III. RESULTS 

Within the T range of the present measurements (2 
K < T < 290 K), a single diamagnetic state (/i + = 
Mu + ) was observed above ~10 K as a signal precessing 
with a frequency uq = 7 M P/27r (comprising the central 
line), where 7 ^ (= 2i r x 135.53 MHz/T) is the muon 
gyromagnetic ratio and B = \B\. The depolarization 
rate for this signal was almost independent of T with 
a Gaussian damping (~ 0.02 /is -1 ); this was unambigu¬ 
ously attributed to random local fields exerted by the 
nuclear magnetic moments of the 47 Ti (natural abun¬ 
dance: 7.4%) and 49 Ti (5.4%) nuclei. Meanwhile, the 
/iSR time spectra exhibited a remarkable change below 
~10 K. A typical example (observed at 5 K, B cx 32.5 
mT) is shown in Fig. [2][a), where a beat pattern due to 
multiple frequency components is clearly observed in the 
signal amplitude. 

As is evident in Fig. [ 2 ][b), one can trace two pairs 
of satellite lines in the fast Fourier-transformed spectra 
(z 7 ±, i = 1 , 2 ), which are situated almost symmetrically 
around the central line (i/q cx 4.41 MHz). A preliminary 
analysis showed that the least-square curve fits of the 
/iSR time spectra, assuming one pair of satellites, did 
not reproduce the data for those obtained with B ori¬ 
ented away from the ( 100 ) direction, yielding a poor % 2 . 
However, fits with two sets of satellites were found to 
yield a satisfactory result with significantly improved x 2 
in the majority of cases. It was also found that the split¬ 
ting {yi+ — Z 7 _) of these satellites remained unchanged 
when the applied field was reduced to 10 mT, indicating 
that B was in the high field-limit range, where the split¬ 
ting was independent of the B value. These observations 
clearly indicate that a significant fraction of implanted 
muons form a paramagnetic Mu state characterized by 
an extremely small HF parameter (almost 10 -4 times 
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FIG. 2: (Color online) (a) Example of /xSR time spectrum 
observed in Ti02 at 5 K under an applied magnetic field of 
32.5 mT, where the field (B) is rotated from (100) axis (a) 
by 0 = 30° around the c axis (9 — 90°). Solid curve shows 
a result of curve fit. (b) Fast Fourier-transform (FFT) of the 
data shown in (a), (c) FFT of the time spectra with B || c at 
5 K. 

smaller than the vacuum value of 4.463 x 10 3 MHz) with 
significant anisotropy in the electronic structure. 

The results of curve fits for the time spectra assum¬ 
ing multiple precession frequency components (z^±) are 
summarized in Fig. [3j for both (a) 0 and (b) 0 depen¬ 
dence. (We assumed that the signal amplitudes of z^ + 
and Vi _ were common for each pair in the curve fits.) 
As can be clearly seen in Fig. [3ja) , Vi± exhibits a strong 
anisotropy with a characteristic angle dependence that is 
symmetric over 0 ~ 45°. Considering the fact that there 
are four inequivalent H (Mu) sites with different princi¬ 
pal axis orientations for anisotropic HF parameters [see 
Fig. D](b)], four sets of satellite pairs discerned by differ¬ 
ent ^-dependence can be expected. More specifically, it 
is predicted that the 0 dependence of the HF parameters 
can be expressed in the form 

Vi±{4>) - "o ± ^A±(4>), (1) 

A±{(j>) = A\ cos 2 $ + A 2 sin 2 <f>, ( 2 ) 

where <f> = 0 db 0 o or 0 d= 0 o =b \, and 0 o is the offset 
angle from the a axis. It must be noted that the dou¬ 
ble sign correspondence between Vi± and the right hand 
side of Eq. O is not unique, because the sign of A_l( 0 ) 
depends on that of Ai and, thus, can change over the 
0 range. While our measurements are insensitive to the 
absolute sign of Ai, they can determine the relative sign 
between A\ and A 2 . The solid curves shown in Fig. [3][a) 
are the result of a simultaneous fit for those curves using 
common values of Ai and 0 o, where A 1 A 2 < 0 is in¬ 
ferred. This strongly suggests that the HF interaction is 
predominantly determined by the magnetic dipolar field 
generated by the off-site electron spin(s). (This scenario 
closely resembles the case of Mu-substituted free radi¬ 
cals often observed in unsaturated organic compounds.) 
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FIG. 3: (Color online) (a) Angle dependence of muon spin 
precession frequency observed at 5 K for (a) azimuthal angle 
(0) and (b) polar angle (6). Filled symbols: Vi± satellite lines, 
open symbols: vq central line. The solid curves are obtained 
via curve fits using general forms of the anisotropic hyperhne 
parameters, while the dashed curves are calculated for a given 
set of parameters (see text). 


Moreover, considering that Eqs. O and m predict eight 
lines when \A±\ 7 ^ \A 2 \, Fig.[3](a) indicates that \Ai \ is al¬ 
most equal to \A 2 \ within the frequency resolution, which 
is limited by the observation time range (~0.1 MHz). 
The obtained values for Ai and 0 are shown in Table H 
Hereafter, we assume that A\ < 0 (see below). 

The 0 dependence is shown in Fig. [3jb) , and is char¬ 
acterized by a narrow splitting of v\± at 6 = 0 [see 
Fig. [21(c)], along with a monotonous increase with increas¬ 
ing 6. Only one pair of satellites has been identified in 
this data set, most likely because of the narrow splitting 
relative to the frequency resolution. It is also notable 
that \v\- — vo | tends to be smaller than |zq + — vq\ for 
some unknown reason. We tentatively attribute this to 
systematic uncertainty in the present measurement, and 
adopt vi+ in order to deduce the HF parameters using a 
form similar to Eq. such that 

A\\ ( 0 ) = As cos 2 0 + A_l( 0 ) sin 2 0, (3) 

where the sign of As is left undetermined. The value of 
As deduced from the curve fit is shown in Table I, where 
|Aj_(0)| has been determined to be 1.314(15) MHz (which 
is in perfect agreement with the value deduced from the 
0 dependence). 
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Mu Ai (MHz) A 2 (MHz) A 3 (MHz) 0 o (deg) 
A±(<f>) -1.29(6) +1.29(6) - 25.5(1.4) 

MO) ~ ~ -0-17(2) 

H -1.276(3)* +1.961(3)* -1.076(3)* 22.1 


TABLE I: Mu hyperfine parameters deduced from curve fits 
of data shown in Fig. [3] (where \Ai\ = |M _2 1 was assumed). 
The A7s for the H-related centei*2i are scaled against the 
gyromagnetic ratio of fi + (indicated by an asterisk, see text). 
The Ai signs for Mu are assumed to be common to those for 
H. 



FIG. 4: (Color online) Hyperfine interactions expected for the 
magnetic dipolar field (Bdip) generated by the spin polariza¬ 
tion of both the Ti and O atoms under an external field (B) 
parallel to the (a) A i, (b) A 2 , and (c) A 3 axes, respectively 
(see text). 


IV. DISCUSSION 


A. Electronic Structure 


It would be intriguing to compare the observed elec¬ 
tronic structure with that of a H-related paramagnetic 
center observed via EPR/ENDOR in a “lightly reduced” 
rutile sample^!. According to the proposed model, the 
center consists of a substitutional Ti 3+ ion adjacent to 
a substitutional [OH] - molecular ion (which is identical 
to the atomic configuration shown in Fig. [TJ apart from 
the Ti valency). The HF parameters for the electron- 
X H nucleus (proton) interaction determined by the EN- 
DOR spectra are reported to be A\ = —0.401(1) MHz, 
A 2 = +0.616(1) MHz, and A 3 = —0.338(1) MHz with 
their principal axes being displaced from the ( 110 ) or 
(110) axes by 22.9° (or by 22 . 1 ° from the a axis). As 
illustrated in Figs. QJc)-(e), the Ai signs are consistent 
with the magnetic dipolar field generated by a d electron 
centered around the Ti site, and their relative magni¬ 
tudes are in reasonable agreement with those expected 
for a point-like dipole at the Ti site, i.e., Ai:A 2 '.As = 
—1:+2: —1. These values are estimated based on the 
diagonal terms of the dipole tensor, A\ oc — |/I e |/r 3 , 
A 2 oc + 2 |/I e |/r 3 , and A 3 oc —[/J e \/r 3 for an effective mag¬ 
netic moment |/J e |, with r being the distance between the 
H and Ti atoms. 

Provided that Mu forms a complex state identical to 
that of H, the corresponding HF parameters can be pre¬ 
dicted by simply scaling these values for H using the fac¬ 
tor 7 M / 7 P = 3.1832 (as shown in Table I, disregarding a 
minor difference of </>o), where 7 P = 27r x 42.5774 MHz/T 
is the proton gyromagnetic ratio. Adopting the reason¬ 
able assumption that the Ai signs are common between 
the Mu and H complexes, comparison of the Ai values 
given in Table I immediately leads us to the conclusion 
that the Mu complex state differs from that of H. 

In the classical limit, the unpaired electron is regarded 
as a point-like magnetic dipole situated at the nearest 
neighboring Ti site, and the HF parameters (per unit 
Bohr magneton, hb) are given by the equation 


7 M 3 cos 2 Ti 
A • = -2^ 


(4) 


where g e is the electron g factor, r (r = |r|) is the co¬ 


ordinate of the Ti 3+ atom with a muon at the origin, 
and Ti is the angle between r and the symmetry axis 
(ti = 73 = 7 r/ 2 , t 2 =0). Assuming that the distance to 
the nearest neighboring Ti 3+ atom r nn = 0.226 nm (i.e., 
an unrelaxed crystal lattice) and that a full 1 ps moment 
(with g e = 2 ) is present, we have A\ n = A™ = — 10.88 
MHz and A£ n = +21.77 MHz, where nn indicates the 
nearest neighbor. These values are far greater than the 
experimental values, which strongly suggests that the 
spatial distribution must be considered for the spin po¬ 
larization, e.g., using 

/l V* [ J 3 / N 3 COS 2 Ti - 1 

A i = 2+ e/XS I d rPs ^ r ’ -2r3-’ 

where r is now the electron coordinate, and p s is 
the spin density corresponding to the difference be¬ 
tween the charge densities of the spin-up and spin- 
down electron(s)22. Comparison of the A™ [obtained 
using Eq. ([5]) for p s (r) = — r nn )] with the ex¬ 

perimental values implies that p s (r) may have a con¬ 
siderable spread to reduce the effective moment size 
Kl ^ VB f d 3 rp s (r)/r 3 . 

Here, we show that both magnitudes and signs of the 
HF parameters can be understood by considering addi¬ 
tional spin polarization of the O atom next to the muon. 
As illustrated in Fig. [4j assuming that spin polarization 
at the Ti site induces small polarization of the ligand O 
antiparallel to the Ti polarization (usually expected for 
the Heisenberg-type exchange interaction), the Ai incor¬ 
porating the contribution of the magnetic dipole at the 
O site are given approximately by 


Ai 

A2 


A 3 





1 2 e 



(6) 


where e is the relative magnitude of the O spin polariza¬ 
tion against that of Ti and ro is the OMu bond length 
(= 0.109 nm). As summarized in Table mi simulated Ai 
for []J e \ = 0.050 ps and e = 0.077 (i.e., e\jl e \ = 0.0039 ps) 
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A 1 (MHz) 

A 2 (MHz) 

N 

S 

00 

l+el « 

Exp. 

-1.29(6) 

+1.29(6) 

-0.17(2) 

- 

Sim. 

-1.29 

+1.09 

-0.17 

0.050/is 0.077 


TABLE II: Comparison of experimental and simulated Mu 
hyperfine parameters assuming antiferromagnetic spin polar¬ 
ization at nn O site, where the experimental values of A\ and 
A 3 were used to evaluate A 2 , |/J C |, and e in the simulation. 


exhibit almost perfect agreement with the experimental 
values. The model also provides a natural explanation for 
the small magnitude of As that the dipole fields from Ti 
and O atoms nearly cancel at the muon site [see Fig. 0(c)], 
suggesting the possibility that the current assumption on 
the As sign may be irrelevant. We note that such residual 
spin polarization of ligand O atoms is usually expected in 
the case of strong d-p hybridization, as has been observed 
for a variety of transition metal oxides^. 

We currently speculate that the difference in the elec¬ 
tronic structure between Mu and H is primarily due to 
the local environment of the specimen. While the present 
Mu complex state was observed in an “as received” crys¬ 
tal, the H-related paramagnetic state was observed after 
a reduction process (annealing at 600° C for 10 min in 
a N 2 atmosphere^). Considering the fact that the re¬ 
duction process also yields O vacancies, it may be of in¬ 
terest to recall the recent investigation based on density- 
functional theory that suggested that H is more stable at 
the O vacancy site, where it adopts a negatively charged 
state to form a Ti 3+ Hy Q complex stated. Thus, it is 
natural to consider the possibility that the paramagnetic 
state observed via EPR/ENDOR might correspond to 
such a state, although the consistency with experimental 
results must be carefully examined. The fact that the 
relative magnitudes of the A{ s for the H-related center 
correspond to the case of e ^ 0 in Eq. © (i .e., A 1 :A 2 :A 3 
= —1:+2: — 1) seems to constitute supporting evidence, as 
it suggests that the relevant H atom has no neighboring 
O atom. 

The fact that the hyperfine interaction is dominated 
by a magnetic dipolar interaction also places a strong 
constraint on the diffusion of in the c channel, because 
the Ai are strongly dependent on the distance to the 
electron (oc 1/r 3 ). The diffusive motion of fi + against a 
stationary d electron at the Ti site would immediately 
lead to the strong damping of satellite signals due to the 
fluctuation of r and hence of A\. 


B. Thermal Property 

The T dependence of the fractional yield (/) for the 
Mu complex state is shown in Fig. [5] The sum of the 
signal amplitude between the Mu complex state and the 
diamagnetic state (not shown) is almost independent of 
T, suggesting that the Mu complex is converted to a dia¬ 
magnetic state above ~10 K. Here, it must be noted that 


/ ~ 0.4 (< 1) does not necessarily indicate two different 
muon sites with different charge states. In non-metallic 
compounds, it is usually expected that the initial yield of 
the paramagnetic state upon muon implantation is pre¬ 
dominantly determined by the density of the short-lived 
free electrons that are produced via radiolysis near the 
end of the muon radiation track. While this density varies 
between compounds, it is independent of T because the 
radiolysis is an athermal process. We assume that the 
muon site in Ti02 is unique (as shown in Fig. [T]d) over 
the entire observed T range, where the initial yield of 
the Mu complex is controlled by the athermal electron 
density, which is independent of T. 

In general, the diamagnetic state can be either posi¬ 
tively or negatively charged. If the electronic energy lev¬ 
els (i?JSy + , E^/°) associated with the Mu complex are sit¬ 
uated near the top of the valence band, a Mu - state can 
be expected at lower T because the Fermi level (Ep) is 
likely to be situated far above the mid-gap in the present 
n-type specimen (i.e., < Ep). In that case, 

an increase in the yield of the Mu complex would be ob¬ 
served with elevating temperature due to the hole capture 
process (Mu - + h + -A Mu); however, this is contrary to 
the actual observed behavior. Thus, the behavior of the 
Mu complex strongly suggests that the process relevant 
to the promotion of the diamagnetic state is electron re¬ 
lease, Mu -A Mu + + e~ (i.e., E^/ + < Ep < E^°), and 
that the Mu complex state can serve as an electron donor. 

Provided that the ionization of the Mu complex state 
is driven by an Arrhenius-type activation process, the 
disappearance of the Mu signal above ^10 K suggests 
that the activation energy (E a ) is of the order of 10 1 
meV. This is in line with certain earlier reports sug¬ 
gesting unidentified shallow level states (e.g, via opti¬ 
cal absorption^ or infrared absorption spectroscopy on 
deuterated rutile^Si). In any case, it must be remembered 
that the interpretation of E a depends on the detailed 
neutral charge state formation process of the Mu com¬ 
plex, and also on the kind of equilibrium state realized 
for the Mu complex formation. At one end of the range, 
E a represents a direct transition from the defect level to 
the bottom of the conduction band, while at the other 
extreme, in equilibrium, it indicates a transition from the 
defect level to the Fermi level (i.e., E a ~ Ep — E^J + ). 
Since the origin of the n-type conductivity in the present 
specimen is unknown, the present value of E a should be 
interpreted as a lower bound for the actual defect level. 


V. SUMMARY AND CONCLUSION 

We have demonstrated that the electronic structure of 
the interstitial Mu center in rutile is characterized by 
extremely small and highly anisotropic hyperfine param¬ 
eters. These parameters are predominantly determined 
by magnetic dipolar interaction with the unpaired Mu 
electron, which is primarily situated at the Ti site. The 
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FIG. 5: (Color online) Temperature dependence of fractional 
yield for muonium complex state deduced from signal ampli¬ 
tude. 

hyperfine parameters are quantitatively explained by a 
model that considers a small residual spin polarization 
of the O atom (which is antiparallel to that of the Ti 
atom), suggesting that the electronic structure should 


be interpreted as being a Ti-O-Mu complex state. The 
extremely small effective moment size of the unpaired 
electron (^0.05 /ib at the Ti site, ^0.0039/i# at the O 
site) as well as the small activation energy required for 
its promotion to the conduction band, implies that the 
Mu complex (and hence the corresponding H state) can 
serve as an electron donor. This strongly suggests that 
H is one of the primary origins of unintentional n-type 
conductivity in rutile TiCU- 
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